Stabilization of the number of Bose-Einstein condensed atoms in evaporative cooling 

via three-body recombination loss 
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The dynamics of evaporative cooling of magnetically trapped ^'^Rb atoms is studied on the basis 
of the quantum kinetic theory of a Bose gas. We carried out the quantitative calculations of the time 
evolution of conventional evaporative cooling where the frequency of the radio-frequency magnetic 
field is swept exponentially. This "exponential-sweep cooling" is known to become inefficient at 
the final stage of the cooling process due to a serious three-body recombination loss. We precisely 
examine how the growth of a Bose-Einstein condensate depends on the experimental parameters of 
evaporative cooling, such as the initial number of trapped atoms, the initial temperature, and the 
bias field of a magnetic trap. It is shown that three-body recombination drastically depletes the 
trapped *^Rb atoms as the system approaches the quantum degenerate region and the number of 
condensed atoms finally becomes insensitive to these experimental parameters. This result indicates 
that the final number of condensed atoms is well stabilized by a large nonlinear three-body loss 
against the fluctuations of experimental conditions in evaporative cooling. 
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I. INTRODUCTION 



Evaporative cooling is the key experimental technique 
in recent successful demonstrations of Bose-Einstein con- 
densation (BEC) with trapped neutral atoms 2, 3j ,4, 
1^ . The cooling mechanism is understood to be that evap- 
oration caused by atomic collisions selectively removes 
the energetic atoms and collisional rethermalization si- 
multaneously lowers the temperature of atoms remaining 
in a trap |^ [3 ■ The utilization of interatomic elastic col- 
lisions makes evaporative cooling highly efficient. Ultra- 
low temperatures of the order of submicro-Kelvin have 
been achieved in BEC experiments by applying evapora- 
tive cooling to laser-cooled atoms. However, in the case 
of real atoms, other undesirable collisions occur, such as 
elastic collisions with background gas and inelastic colli- 
sions due to dipolar relaxation and three-body recombi- 
nation [1,0|- These collisions result in the loss of trapped 
atoms and seriously reduce the effectiveness of evapo- 
rative cooling. The efficiency of evaporative cooling is 
therefore determined by the competition between evapo- 
ration and the loss due to the undesirable collisions. 

Several theoretical analyses based on a classical gas 
model have addressed the complicated dynamics of evap- 
orative cooling 6,i 8,1,,10, 11, 12, 13, 14, 15, 16, 17]. 
The results have been applied to the optimization of ex- 
perimental conditions u3, Il7l| and have provided useful 
guidelines for reaching the BEC transition point experi- 
mentally. However, the classical theories fail in the quan- 
tum degenerate region at low temperatures. Bose statis- 
tics of atoms strongly modifies the dynamics of evap- 
orative cooling there: interatomic elastic collisions are 
enhanced by the stimulated scattering process, and the 
event rates of inelastic collisions due to dipolar relax- 
ation and three-body recombination are reduced by the 
coherence of identical condensed atoms d&j 

Re- 
cent theoretical works HH HHH H El lH have proved 
that the quantum kinetic theory taking into account the 



bosonic nature of atoms is a powerful tool to gain a com- 
prehensive understanding of evaporative cooling in BEC 
experiments. 

In our previous paper |27| . we studied the optimiza- 
tion of evaporative cooling in a magnetically trapped 
^^Rb system on the basis of the quantum kinetic the- 
ory. We demonstrated that the acceleration of evapo- 
rative cooling around the BEC transition point is very 
effective against the loss of trapped atoms due to three- 
body recombination. The number of condensed atoms is 
largely enhanced by simply optimizing a sweeping func- 
tion of the frequency of radio-frequency (rf) magnetic 
field. On the other hand, this optimized cooling is deter- 
mined for a fixed certain experimental condition and its 
cooling trajectory is sensitive to the experimental param- 
eters of evaporative cooling, such as the initial number 
of trapped atoms and the initial temperature. The opti- 
mized evaporative cooling therefore might show the fluc- 
tuations of the final number of condensed atoms due to 
some variations of these experimental parameters inher- 
ent in the real experiments. In current BEC experiments, 
the number of condensed atoms is usually measured by 
the destructive time-of-flight method. The stable pro- 
duction of condensates with small number-fluctuations is 
thought to be another essential property of evaporative 
cooling. 

In this paper, we quantitatively analyze the dynam- 
ics of the conventional evaporative cooling where rf- 
frequency is swept exponentially. This "exponential- 
sweep cooling" is known to become inefficient at the 
last stage of the cooling process due to three-body re- 
combination loss. We precisely examine the dependence 
of condensate growth on the experimental parameters, 
such as the initial number of trapped atoms, the initial 
temperature, and the bias field of magnetic trap. We 
demonstrate that three-body loss drastically depletes the 
trapped ^^Rb atoms around the BEC transition point 
and that the number of condensed atoms finally becomes 
insensitive to these experimental parameters. The final 
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number of condensed atoms is well stabilized by a large 
nonlinear three-body loss against the fluctuations of ex- 
perimental conditions. 



collisions. The change rates of the total number of atoms, 
N, and the total internal energy, E, are thus calculated 
as the sum of the contributions of all processes |23| '■ 



II. QUANTUM KINETIC THEORY OF 
EVAPORATIVE COOLING 

In this section, we briefly mention our theoretical 
framework (See Ref. H^l). During evaporative cooling, 
the magnetic trapping potential U{r) is truncated at the 
energy, which is determined by the frequency of the ap- 
plied rf-magnetic field. The thermalized distribution of 
the noncondensed atoms in such a truncated potential is 
well described by the truncated Bose-Einstein distribu- 
tion function such that 



f{r,p) = 



1 



exp(ep/fcBr)/e(r) - 1 



Q{A{r)-ep), (1) 



where = p^/2rn is the kinetic energy of atoms with 
momentum p and mass m, T is the temperature, and 
0(a;) is the step function. We add tilde to the nota- 
tions of the quantities obtained through this truncated 
distribution function. The function ^(r) represents the 
local fugacity and includes the mean-field potential en- 
ergy such that ^(r) = exp{[^ — C/(r) — 2vh{r)]/kBT}, 
where /i is the chemical potential, v — Airah^ the interac- 
tion strength of atoms in proportional to the s-wave scat- 
tering length a, and n(r) the number density of atoms. 
The step function in Eq. eliminates the momentum 
state whose kinetic energy exceeds the effective poten- 
tial depth A{t) = et — U{r) — 2vh{r), where Ct is the 
truncation energy of a magnetic trapping potential. The 
number density of atoms n{r) and the internal energy 
density e(r) are evaluated through the truncated distri- 
bution function in a self-consistent way: 



h{r) 



/(r,p) p^dp, 



in 

^ J '^P p'^dp + v[ii{r)f 

+C/(r)n(r). 



(2) 



(3) 



Thus both h{r) and e(r) are the complicated functions 
of T, /i, ei, U{r), and a. The spatial integrations 
of these density functions give us the total number of 
atoms, N — J n{r)dr, and the total internal energy, 
E = J e{r)dr. After the BEC transition, the number 
density in the condensed region is given by the sum of a 
condensate no{r) and saturated noncondensed atoms n„ 
such that n{r) = no{r) + hn- The condensate no{r) is ap- 
proximated by the Thomas-Fermi distribution and fin is 
evaluated from the truncated Bose-Einstein distribution 
function on the condition that the local fugacity becomes 
unity, i.e., |(r) = 1. 

During the cooling process, trapped atoms are removed 
from the potential by both evaporation and undesirable 



dN 
It 



dE 
'dt 



j ncv(r) dr-Y, GsJ K.(r)[n(r)]^ dr 



dN 



(4) 
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eev(r) dr -Y,Gs if.(r)g(r)[n(r)]^-i dr 

s=l 

(5) 




The ncv(r) and ecv(r) denote the evaporation rates of 
density functions derived from a general collision inte- 
gral of a Bose gas system. In Eqs. Q and we 
adopt the opposite sign of the notations of these evapo- 
ration rates, which were defined in our previous paper 
[23| . Parameters Gi, G2, and G3 are the decay rate 
constants of trapped atoms due to the background gas 
collisions, dipolar relaxation, and three-body recombi- 
nation, respectively. The function Ks(r) represents the 
correlation function, which describes the s-th order co- 
herence of trapped atoms, and we use the expressions 
of Ks for an ideal Bose gas system 0, |23|- The 
terms proportional to the change rate of truncation en- 
ergy, it = det/dt, give the contribution of extra atoms 
that "spill over" when et decreases continuously in forced 
evaporative cooling [l6ll27|. 

The time-evolution was calculated numerically by as- 
suming that the system always stays in quasi-thermal 
equilibrium state and is described by the truncated Bose- 
Einstein distribution function. This quick rethermal- 
ization approximation is appropriate for the slow evap- 
orative cooling normally adopted in BEC experiments 

MM. 



III. NUMERICAL RESULTS AND DISCUSSINS 

In this section, we show the numerical results of the 
time-evolution calculations of evaporative cooling corre- 
sponding to the several cases of experimental conditions. 
To carry out the quantitative calculations, we employed 
the relevant parameters corresponding to the experimen- 
tal setup at Gakushuin University The cloverleaf 
magnetic trap characterized by bias field Bq, radial gra- 
dient B , and axial curvature B was modeled by a mag- 
netic field of the loffe-Pritchard type such that B{r, z) = 
^(ar)2 -f {f3z^+BoY, where a = yjB"^ - B"Bq/2 and 
fi = B /2. We used the typical values adopted exper- 
imentally, such as i?o = 1-6 G, B = 169 G/cm, and 
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B = 148 G/cm^. The corresponding trapping frequen- 
cies are = 27r x 170 Hz in the radial direction and 
ujz — 2tt X 15.5 Hz in the axial direction. 

The frequency of the rf-magnetic field, I'lf , is assumed 
to be swept exponentially from 35 to 1.13 MHz in 30 s 
such that J^rf(i) = i^a + — Va) 6^*/'^ with Va ~ 1.13 
MHz, Vb = 35 MHz, and r = 3.69 s. This exponential- 
sweep cooling is commonly adopted in EEC experiments 
[20j . and we employed this functional form in the fol- 
lowing all numerical calculations. Furthermore, we set 
the collisional parameters of ^^Rb atoms in the trapped 
state of F = 2,mF = +2 as follows: s-wave scattering 
length a = 5.8 nm, loss rate constants due to collisions 
with background gas Gi = 0.01 s^^, dipolar relaxation 
G2 = 1.0 X 10~^^ cm^s"""^ [23, and three-body relaxation 
G3 = 1.1 X 10-28 cm^s-iTlg. 



A. Dependence on the initial number of trapped 
atoms 

Considering that the initial temperature is strongly 
restricted by the limit of the laser cooling done before 
evaporative cooling '3^, we first carried out the time- 
evolution calculations for several initial numbers of atoms 
with the initial temperature fixed. Under the typical ini- 
tial temperature of Ti = 380 /j,K |2^ , we varied the initial 
total number of atoms from Ni = 1.5 x 10^ to 3 x 10^ and 
investigated how the final number of condensed atoms 
depends on these values. 

Figures ^a) and 1(b) respectively show the time- 
evolution of the total number of atoms N and tempera- 
ture T for several initial conditions of evaporative cool- 
ing. We can see that both N and T decrease almost 
exponentially over time. Furthermore, in Fig.^a), four 
of the curves from the solid line to the dash-dotted line 
corresponding to A^^ = 3 x 10^ ~ 3 x 10^ merge into 
one, which means ten times more initial atoms are com- 
pletely lost during the 30-s cooling. The time-evolution 
of temperature in Fig. Q^b) also shows the same feature. 
The final results of exponential-sweep evaporative cooling 
become insensitive to the initial number of atoms when 
cooling starts from more than 3 x 10* atoms. 

Figure[21shows the dependence of the time evolution of 
the number of condensed atoms, A^Oi on the initial num- 
ber of trapped atoms Ni. For JV; = 3 x 10* ~ 3 x 10^, 
we finally obtain the same A^o value of 8.3 x 10^ with 
the same condensate fraction of Nq/N = 51% after 30-s 
evaporative cooling, even though the BEC transition oc- 
curs at a different time. This implies that, when we apply 
exponential-sweep evaporative cooling, even a one-order 
change of the initial total number of trapped atoms does 
not affect the final number of condensed atoms. Since 
A^j = (3 ~ 5) X 10* is normally achieved in current 
BEC experiments with *^Rb, it is concluded that the 
final number of condensed atoms can be well stabilized 
against the fluctuations of the initial number of atoms in 
the experiments. 
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FIG. 1: Time evolution of (a) total number of atoms and 
(b) temperature T for several initial conditions of evaporative 
cooling: the initial number of trapped atoms Ni is varied from 
1.5 X 10* to 3 X 10® under the fixed initial temperature of 
Ti = 380 /iK. In both figures, the curves correspond to cases 

with Ni = 1.5 X 10* ( ), 2 x 10* ( ), 3 x 10* ( 

), 5 X 10* ( ), 1 X lO'' ( ), and 3 x 10® ( ). 
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FIG. 2: Time evolution of the number of condensed atoms 
Nq for the several initial numbers of trapped atoms Ni. We 
use the same lines as in Fig. in terms of the corresponding 
Ni values. The solid line (A^i = 3 x 10®) and the dotted 
line (Ni = 1 x 10®) almost coincide with each other. For 
Ni — 1.5 X 10*, the system can not reach the BEC transition 
point and condensed atoms are not produced. 

Here we show the results of the cooling trajectory 
through phase space and provide clear information about 
the efficiencies of present exponential-sweep cooling. In 
a nonuniform gas system, the phase-space density is 
evaluated at the peak position as p = n(0)A'^, where 
X — hj \/2-KmkBT is the thermal de Broglie wavelength 
with atomic mass m. In Fig. |31 we plot the cooling 
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FIG. 3: Cooling trajectories through phase space for several 
initial conditions of evaporative cooling: the initial number of 
trapped atoms Ni is varied under the fixed initial temperature 
of Ti = 380 fiK. The phase-space density p is normalized by 
the critical value pc = 2.612. The BEC transition occurs when 
each trajectory passes the horizontal line where p/ pc = 1. We 
use the same lines as in Fig.Qin terms of the corresponding Ni 
values. The arrow shows the direction of time. The symbols 
on the solid line (circle), dotted line (triangle), and dashed 
fine (square) indicate the points where the evaporative cooling 
becomes significantly inefficient. 



trajectories of the present evaporative cooling in the iV- 
p plane. At the early stage of coohng, all trajectories 
are well separated according to the initial conditions 
and evaporative cooling works efficiently. However, as 
the system approaches the quantum degenerate region 
where the phase-space density exceeds the critical value 
of Pc — 2.612, three of the trajectories from the solid 
line to the dashed line bend towards horizontal around 
the respective turning points indicated by the symbols, 
showing the cooling efficiency drastically decreases there. 
As expected from the results in Fig. ^ three of the tra- 
jectories merge and finally coincide with the dash-dotted 
line corresponding to iVi = 3 x 10^. On the other hand, 
when Ni < 2x 10^, the final results after 30-s evaporative 
cooling strongly depend on the initial number of atoms. 
The system can not reach the BEC transition point if the 
cooling is started from Ni = 1.5 x 10^. 

Next, we show the loss rates of trapped atoms during 
evaporative cooling and discuss the origin of the ineffi- 
ciency of cooling. As described in Eq. I^J, the change 
rates of the total number of atoms decompose into five 
contributions: the evaporation rate N^v^ spill rate A'spiu, 
loss rates due to back-ground gas collision Ni, two-body 
dipolar relaxation iV2, and three-body recombination N^. 
In Fig. 21 we plot the time evolution of these rates for 
iV, = 3 X 10^ 1 X 10^ and 5 x 10®. The results for iVi are 
not shown here since this rate is given by GiN and evalu- 
ated straightforwardly from Fig.^a). The symbols in the 
figures correspond to the same points plotted on the cool- 
ing trajectories in Fig.|21and indicate the time where the 
efficiency of evaporative cooling is drastically reduced. 
We see that both Ncv and Nspm [U decrease exponen- 



tially over time, while N2 and A^3 have peak structures 
around the symbol on each curve. The peak value of N^ 
is about five times larger than that of N2 and three-body 
recombination is the dominant loss mechanism for ®^Rb 
atoms. By comparing Fig.^fa) with^fd), it is also found 
that the symbol on each curve indicates the time when 
the three-body loss rate exceeds the evaporation rate. At 
the final stage of the present exponential-sweep cooling, 
we understand that the cooling speed becomes so low 
that the three-body loss crucially reduces the cooling ef- 
ficiency. A stabilization of the final number of condensed 
atoms is therefore realized by the combination of slow 
cooling and large three-body loss. 

From another viewpoint, the present result also sug- 
gests that the exponential-sweep of rf-frequency is se- 
riously inefficient when we try to make a larger conden- 
sate. We have thought so far that the improvement of the 
initial conditions results in more condensed atoms even 
in exponential-sweep evaporative cooling. However, it is 
clear in Fig. 13 that a ten- fold increase of the initial num- 
ber of trapped atoms does not improve the achievable 
number of condensed atoms. In order to obtain a larger 
number of condensed atoms, we have to carry out the 
evaporative cooling quickly at the final stage, where the 
system approaches the quantum degenerate region with 
high density, as has been demonstrated quite recently by 
optimizing the cooling trajectories |15Hl'<if2^ . 

B. Dependence on the initial temperature 

In BEC experiments, evaporative cooling is normally 
employed after laser cooling ^ .3Qj . The atoms are trans- 
fered from a magneto-optical trap (MOT) in laser cooling 
into a magnetic trap in evaporative cooling. The mis- 
alignment between these two traps easily causes the fluc- 
tuations of initial temperature for succeeding evaporative 
cooling process. Thus, we next carried out the time- 
evolution calculations of evaporative cooling for several 
initial temperatures with the initial number of trapped 
atoms fixed. We assumed that 5 x 10^ atoms are ini- 
tially loaded into the magnetic trap before evaporative 
cooling (i.e., 5 x 10^ atoms exist in the magnetic trap 
without truncation). The initial temperature was then 
varied from Ti — 280 p,K to 630 p,K, and we investigated 
how this variance of initial temperature, as well as that of 
the initial number of atoms, is stabilized by a three-body 
loss. 

Figures |3Ia) and \^h) respectively show the time- 
evolution of the total number of atoms N and the temper- 
ature T for several initial conditions of evaporative cool- 
ing. Both N and T decrease almost exponentially over 
time. In Fig.[S^a), the small difference of iV at i = orig- 
inates from the difference of truncation-parameter value 
77 = et/{kBTi) for each initial temperature with the same 
truncation energy ct- We can see that N curves and T 
curves for the initial temperatures, i.e., from 280 /iK to 
580 /iK, finally merge into one. The variance of initial 
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FIG. 4: Time evolution of (a) evaporation rate A'cv, (b) spill rate A'spin, (c) loss rate due to two-body dipolar relaxation N2, 
and (d) loss rate due to three-body recombination Ng. We use the same lines as in Fig. 0in terms of the corresponding Ni 

values, and show the results for the cases of A^i = 5 x 10* ( ), 1 x 10^ ( ), and 3 x 10'' ( ). The symbols in the 

figures correspond to the same points plotted on the cooling trajectories in Fig. |H1 



temperatures finally vanishes after 30-s exponential evap- 
orative cooling as similarly in Fig. ^ 

In Fig. El we plot the growth of condensate at the final 
stage of cooling. For Ti — 280 ~ 580 /zK, the final num- 
ber of condensed atoms at i = 30 s takes almost the same 
value, A^o — 8.3 x 10"*, even though the BEC transition 
occurs at different times according to the initial tempera- 
tures. We will obtain the same number of condensate for 
a wide range of initial temperatures. When Ti exceeds 
580 /iK, on the other hand, the final value of Nq de- 
creases with increasing Ti. It is concluded that the num- 
ber of condensed atoms is sufficiently stabilized against a 
large variance of initial temperatures too (300 /^K in the 
present calculation). We also verified that this stabiliza- 
tion is caused by a combination of large three-body loss 
and slow evaporation at the final stage of cooling process, 
which is exactly the same mechanism as discussed in the 
previous subsection. 

C. Dependence on the bias fleld of a magnetic trap 

In BEC experiments, the magnetic trap is usually 
driven with a large current of the order of 100 A. The 
current noise causes the fluctuations and drifts of the 
bias field Bq in a magnetic trap, which is the important 
limitation to the stable production of condensates. For 
example, Ref. reported that the dirift of Bq should 



be suppressed to less than 3 niG to obtain a good accu- 
racy for the number of condensed atoms produced exper- 
imentally. Here we study the dependence of condensate 
growth on the bias field of a magnetic trap. The time- 
evolution calculations of 30-s evaporative cooling were 
carried out for several values of Bq ranging from 1.56 to 
1.614 G, under the initial conditions of iVi = 5 x 10^ 
and Ti = 380 /iK. Since we used the same exponential 
sweeping- function of rf-field frequency frf (t), the time de- 
pendence of a truncation energy, et(t), was equivalent in 
these calculations. Therefore, at any time during evapo- 
rative cooling, a larger bias field corresponds to a smaller 
depth of a truncated trapping potential. 

Figure [7| shows the time evolution of the number of 
condensed atoms Nq for the several values of bias field. 
The condensate growth curves strongly depend on Bq 
values as has been expected. With larger bias field, 
the BEC transition occurs earlier and the curve shows 
a peak structure. A condensate gradually decays due to 
three-body loss since slow evaporation at the last stage of 
cooling can not preserve the large number of condensed 
atoms. On the other hand, with smaller bias field, 30-s 
evaporative cooling is completed before the condensate 
grows sufficiently. In Fig.|Sl we plot the bias field depen- 
dence of the final number of condensed atoms Ng-p. It 
is found that TVof is the convex function of Bq and the 
stationary point exists around Bq — 1.585 G. When Bq is 
set at this value, we can expect that the final number of 
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FIG. 5: Time evolution of (a) total number of atoms N and 
(b) temperature T for several initial conditions of evaporative 
cooling: the initial temperature Ti is varied under a fixed 
initial number of atoms A'^i ~ 5 x 10*. In both figures, the 

curves correspond to the cases with Ti — 630 fiK ( ), 

580 fiK ( ), 480 iiK{ ), 380 /iK ( ), and 280 /iK 

( )■ 
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FIG. 6: Time evolution of the number of condensed atoms 
A^o for several initial temperarures of evaporative cooling. We 
use the same lines as in Fig.j^in terms of the corresponding 
Ti values. 



condensed atoms is stabilized against fluctuations of bias 
field within about 10 mG. Note that this stabilization is 
also achieved by a large three-body loss. 



IV. CONCLUSION 

On the basis of the quantum kinetic theory of a Bose 
jas, we have theoretically studied the dynamics of evap- 
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FIG. 7: Time evolution of the number of condensed atoms 
A^o for several values of bias field. The curves correspond to 

cases with Bo = 1.61 G ( ) 1.60 G ( ), 1.59 G ( ), 

1.58 G ( ), 1.57 G ( ), and 1.56 G ( ). 
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FIG. 8: Bias field dependence of the final number of con- 
densed atoms, A^oF. Around the stationary point of Bo = 
1.585 G, we can see that A^of is sufficiently stabilized against 
fluctuations of the Bo value within about 10 mG . 



orative cooling with magnetically trapped ^^Rb atoms. 
The time-evolutional calculations were carried out quan- 
titatively for the conventional exponential-sweep cooling, 
which is known to be inefficient against a three-body re- 
combination loss. We precisely examined the dependence 
of a condensate growth on the experimental parameters, 
such as the initial number of atoms, the initial tempera- 
ture, and the bias field of a magnetic trap. It was shown 
that three-body recombination drastically depletes the 
trapped atoms as the system approaches the quantum 
degenerate region and the number of condensed atoms 
finally becomes insensitive to these experimental param- 
eters. The final number of condensed atoms is well stabi- 
lized by a large nonlinear three-body loss against fluctua- 
tions of experimental conditions. We can choose the man- 
ner of evaporative cooling according to the aim of BEG 
experiments; a conventional exponential-sweep cooling if 
we want a stable production of condensates or an opti- 
mized one with the rapid sweep of rf-frequency at the 
final stage of cooling procedure if we want to maximize 
the number of condensed atoms 1271. 
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